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Introduction 

In  [1]  M.J.  Lighthill  proposed  a  new  perturbation  tech- 
nique for  the  approximate  solution  of  certain  types  of  non- 
linear differential  equations  near  their  singular  points. 
This  method  has  since  found  numerous  useful  applications 
in  hydrodynamical  problems,  an  excellent  account  of  which 
can  be  found  in  a  report  by  H.S.  Tsien  [2],   The  procediores 
are  justified,  in  the  existing  literature,  by  very  convinc- 
ing plausibility  considerations  but  mathematically  complete 
proofs  seem  to  be  lacking. 

The  purpose  of  this  note  is  to  give  such  a  proof  in  the 
simple  but  typical  case  of  the  differential  equation 


(X  +  au)^  +  q(x)u  -  r(x)  =  0,  (1.1) 


where  a  is  a  small  positive  parameter.   This  is  the  example 
used  also  in  [1]  for  an  introduction  to  the  technique.   In 
fact,  our  proof  is  essentially  a  completion  of  the  order- 
of -magnitude  arguments  employed  in  that  article.   We  shall 
therefore  be  somewhat  brief  here  in  our  explanation  of  the 
formal  aspects  of  the  procedure. 
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The  problem  is  to  represent  the  solution  of  (1.1) 
through  a  point 

X  =  a,  u  =  b  (1»2 ) 

in  a  form  that  remains  valid  near  x  =  0,   a  =  0.  This  is 
achieved  by  referring  x  and  u  to  a  suitably  chosen  para- 
meter 4,  as  will  be  explained  in  the  next  section. 

2 •   Construction  of  the  solution . 

Ve   make  the  following  assumptions  concerning  the  dif- 
ferential equation  (1.1)  and  the  prescribed  initial  values 
(1.2)  « 

(i)     a>0,     a>0,      b>0.  (2.1) 

(ii)    The  functions  r(x)  ,  q(x)  are  holomorphic  in  |xl<a  . 
Without  loss  of  generality  we  may  set 

a  =  1   ,  {2,2) 

since   this   can  be  brought  about   by   the    transformation 
X  =  x""'a    ,    a  =   a''""a    . 
VJe   set,    formally. 


2 
X  =   4  +  x,a  +  x„a     +    ...  (2.3) 


2 

u  =  u^  +  u^a  +  Uptt     +    ...      ,         (2,14.) 
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where  the  u  and  x^  are  fxinctlons  of  £  to  be  determined  and 
insert  this  into  the  equation 

(x  +  au)u  +  (q(x)u  -  r(x))x  =  0  ,  (2.5) 


equivalent  with  (1*1)  •   If  r(x)  and  q(x)  are  developed  in 
powers  of  x  -  ^  about  the  point  4  ,  setting 


r(x)  =  r(4)  +  r'(^)(x  -  ^)  +  |  r"(4)(x  -  ^)^   -f 


etc«,  the  left  member  of  (2,5)  becomes,  after  collection  of 
like  powers  of  a,  a  series  of  the  form 

2 
Ij^   +  L-j_a  +  L2a  +  ...  ,  (2.6) 

whose  coefficients  are  fairly  complicated  expressions  in 
the  X.  ,  u.  ,  X.  ,  u.   and  ^  •  (Ue  denote  differentiation 

J         J         J         «i 

with  respect  to  ^  by  dashes.)  More  precisely,  we  find 
(writing  q,r  for  q(^),  r(£)). 


Lq  E  CU^  +  qu^  -  r  ,  (2.7) 

L,5  4^1  +qu  +  (qu^-r)x'  +  (u'  +q'u^-  r'  )x^  -  R^,   s>0    (2.8) 


where 
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the  function  h  being  defined  as  follows:   In  the  expression 
h(x,u,x',^)  =  |-[q(x)  -q-  (x-4)q']u  +  [r(x) -r  -  (x-4)r']|x  (2. 


10) 


replace  x  and  u  by  the  series  (2,3),  {2,k),    expand  about  x  =  ^, 

and  collect  like  powers  of  a  ,   Then  the  coefficient  of  a^ 

is  called  h  .   The  quantity  h  is  therefore  a  polynomial 
s  s 


in  the  x.  ,  x.  ,  u,  whose  terms  are  of  the 


form 


ACOx^-'  x^""      ...  x^^"  x^*"  u^   ,  (2.11) 


with 


a,  +  2ap  +  .  .  .  +  |ji,a  +  vP  +  py  =  s   ,  (2.12) 


P  =  0,  1    ,      Y  =  0,  1 


Since  the  expansion  of  (2,10)  in  powers  of  x  -  ^  contains 

no  constant  or  linear  terms,  we  know  also  that  "^   a.  >  2  . 
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It  follows  that,  if  a  ,  p,  or  y  are  positive,  then 


H<s-1,    v<s-l,    p<s-l. 

In  particular,  h,  is  identically  zero. 

In  order  to  satisfy  the  equation  (1.1)  the  functions 

X  (f )  ,  u  (^)  must  now  be  successively  determined  in  such  a 
s       s 

way  that  the  L   ,  as  defined  in  (H.7)  and  (2,8),  vanish.   At 
the  same  time  we  must  impose  boundary  conditions  on  these 
functions  such  that  resulting  series  {2,3),    (2,1+)  assume 
the  Values  x  =  1  ,  u  =  b  for  some  ^  =   E,-,    *      Now,  each  equa- 
tion L  =  0  ,  s  >  0  ,  can  be  replaced,  in  infinitely  many 
ways,  by  two  differential  equations  of  the  form 


(qu^  -  r)x'  +  (u'  +  q'u„  -  r')x„  =  Rf-*-^         (2.13) 
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C  u^  +  qUg  =  r(2)  (2.11^) 


the   only  condition  being  that 


R^^^   4  r(2)      =  R^      .  (2.15) 
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This  is  a  genuinely  recursive  system,  since  R   depends  only 
on  the  X.  ,  u.  with  j  <  s-1  • 

J         J 

The  radius  of  convergence  of  the  resulting  series  («f,3), 

{2,1^.)   depends  on  E,   and  on  the  manner  in  which  each  R  is 

split  into  R^^  and  R^  ^   •   In  the  next  section  it  will  be 
s        s 

shown  that  an  acceptable  method,  from  this  viewpoint,  consists 

(P ) 

in  taking  R^  '  =  0  for  all  s.   But  there  are  other  possibi- 

lities,  which  may  offer  computational  advantages. 


3«   The  Linear  Recursion  Equations. 

In  order  to  calculate  u  from  the  equation  L  =  0,  (cf , 

o  ^        o    ' 

formula  (2.7))  we  observe  first  that  for  q  =  q(0)  4  0  the 
corresponding  homogeneous  equation  is  solved  by 


-  rq(C)€"^d^     -% 
e  •*  =4    P(4) 


'o 


where  p{£,)    is  holomorphic  and  different  from  zero  in 
Ul  <  1  .   Hence, 

%(^)  =  ^o^  %(^)  +  ^  ""^  P(0  (   p-l(t)t^^"^(t)dt.   (3.1) 

^^1 

The  constant  c^  depends  on  the  choice  of  u  (£, )  •   For  the 
o   -^  o  ^1 

sake  of  simplicity  we  set 
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^1  =  1 

u^d)   =  b   ,      Ug(l)  =   0;      s   >  0  (3.2) 

Xg(l)   =   0   ,  s   >  0  (3.3) 


and  accordingly 

c^     =  b/p(l)    . 
Then    the  resulting  functions  x(^),    u(4)  will  satisfy   the 
boundary  conditions   x(l)   =  1,    u(l,)   =  b    . 

It   is    seen  from   (3 •I)    ttiat  u     has   the   form 


u^(^)   =   4     \iU)   +  %2^£.)  (3*U) 


with  u  , (4),  ^qp^^)   holomorphic  in  \e\    <  1 
In  this  section  we  shall  assume  that 


qQ  >  0  .  {3*S) 


The  case  q^  <  0  will  be  discussed  in  section  6. 
^o 

Next  we  calculate  x,  .   Inspection  of  formula  (2.9) 

shows  that  R,  reduces  to  u  u   .   Following  Lighthill  we 
1  o  O  o    o        — 

require  that  Ri  '   contain  the  part  of  u  u  that  involves 
negative  powers  of  g  higher  than  g    ,   '-i-hen  we  have 
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-,,      -2q  -1 
H       =   ^       ^ll(^)    '  (3.6) 

where  r,-,(4)  is  a  polynomial  in  4    with  coefficients  that 
are  holomorphic  in  4  for  |^|  <  1  .   Inserting  (3 •6)  and  the 
expression  for  u  and  u  into  (2,13)  that  equation  is  seen 
to  be  of  the  form 


^e{i)x^  +   f(C)x-L  =   4    r^^(4)  (3.7) 

where  e(4)  and  f(4)  are  linear  functions  of  4    with  holo- 
morphic coefficients.   Furthermore, 

ft      1      » 

m  =  ^ ' .u, -^  -  =  i.^a^o(,)) 


=  (^  -  ^^(1  -^  0(^^)  =  -  1  +  0(4)  +  0(4^°) 
^  o 


At   this  point  our  procedure  breaks  down  unless  we  impose 
a  condition  not  mentioned  in  [1]  or  [2],  viz. 

q(?)  %{£>)   -  r(4)  =1=0,      inO<4ll.  (3.8) 

Then  e(4)  does  not  vanish  in  this  interval.   Hence, 
exp  I   -    (  tM,   d4  \   =   471(4)   , 
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where  ti{£,)    Is,  in  0  <  ^  <  1,  a  holomorphlc  function  of  the 
two  variables  4  and  ^    without  zeros. 

The  solution  of    (3.7)   with  x^(l)   =   0  is 

x^(4)=-47i(4)  ^    u"^(t)t     °       r^^(t)dt  =  c     °x^'(£).  (3.9) 


-1/^N^   ^o        _      ,^vj^      ^  ^o   -:s- 

71  -  - 
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The  function  xi!{E,)    is  holomorphic  in  ^  and  ^  °  , 

For  a  fuller  description  of  the  stepwise  calculation  of 

the  X  (4)  and  u  (4)  we  refer  to  [1]  and  [2]  .   At  each  step 
s         s 

(2  ) 
R  must  he  split  in  such  a  way  that  R^  '  is  of  the  order 

-(s+l)q  n) 

0(4        )   as  4  -^  0.   Any  such  choice  of  R*-  '   and 

(2  ) 
R^  '   will  be  called  admissible.   Then 
s 

(1)  -(s+l)q    -1 

R^^^     =0(4  °      )     ,  (3.10) 

and 

-s^o     ....  -(s+l)q        ,, 

^3(4)  =4     ^  Xg(4)   ,   u^(4)  =  4  ^'sU)  (3.11) 

q 

the  fiinctions  x'"'(4)    ,    u"(4)     being  holomorphlc    in   4  and  4        . 

s       s 

These  factg  suggest  the  conjecture  that  the  resulting 
series  for  x(4)  and  u(4)  have  radii  of  convergence  of  the 
order  0(4    )  .   However,  in  view  of  the  degree  of  arbitrari- 
ness in  the  decomposition  of  R   into  R,^  '   and  R^  '  one  can 

S  X  o 

expect  to  prove  such  a  convergence  statement  only  on  the 


^  .'     .  :\    tirtl'A    ,  \  ^c  i     to    L-i. 


>       M. 


■:..:•'*   . '.:    VOt-jQ  ■  ■'  yi:..,:. 


::>;.£  f!    .)/,■!■; -^    .".'. 


1       j,.-y  J  -.J  * 


(.:    l;n^      .   .i  '  ,.   ■     •-•>■ 


uT     \-,t>f    i-    iijjj? 


I  »/ri  » 


n::^f    ^.;I, 


,  >  .1.  t'^  ! 


{■:    * 


5    '   ■'•: 


•>-       \  « 


'  B 


..■<>; 


-i         i       '1'  Z'; 


f  ^ 


ons 


.}    :   .    V.'  .(<■    J,.    ;-  '^;v./  ? 


-  10  - 

basis  of  some  specific  rule  for  this  decomposition.   If  we 

wish  to  keep  the  choice  of  R^^  free,  within  the  limitations 
^  s 

imposed  above,  we  must  be  content  with  a  iitatement  of  asymp- 
totic character. 

The  method  to  be  used  here  is  a  combination  of  these  two 
possibilities.  Vie  shall  admit  any  admissible  choice  of  R^  ' 
for  s  <  m  ,  m  arbitrary.   For  s  >  m  we  shall  set,  however, 

rI-"-^   =  R.   ,    rI^^   =0   ,   s  >  m  .     (3.12) 
s         s        s 

Hence,  the  series  (2,[^)  terminates  at  the  term  in  a   ,  thanks 
to  our  choice  of  initial  conditions, 

AS  a  preliminary  to  the  convergence  proof  we  derive  two 
simple  inequalities  concerning  the  solutions  of  the  diffe- 
rential equation  problem 

?e(4)v'  +  f(4)v  =  <|)(4)   ,   v(l)   =   0  ,   (3.13) 

where  e(^)  and  f(,^)   have  the  same  meaning  as  in  (3*7),  and 
^{^)    is  uniformly  bounded  in  every  interval  0  <  B,     <   E,  <  1    . 

The  solution  v(f )  is  given  by  the  formula  (3*9)  with  rTT(t) 

q 
replaced  by  -  (j)(t)t  °  .   If 
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0  <  m  <  |7t(4)|  <  M  ,   for  0  <  4  1  1   »  (3.1U) 

we  find 

\v(^i<  £,   ^  sup   |<t)(t)|(4"'^-  1)  <csup   |<t)(t)|,  0  <  4  1  1, 
€5t<l  4<t<l 

M 
provided  c  >  —  .   If  this  result  is  combined  with  (3.13) 

we  have,  making  also  use  of  assumption  (3.8),  that 


|v'(4)|  <   cC'^   sup    \^{t)\      ,      0  <  ?  <  1   , 
C<t<l 


if  c  is  chosen  sufficiently  large. 

V'e  apply  these  inequalities  to  equation  (2,13)  which  is 
of  the  form 

obtaining  the  inequalities 
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|x!(4)l  <  c  4  °   sup   1r^^)'  (3.16) 
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i| •      The  Convergence   Proof  for   q(0)   >  0    . 

Except  for  a  few  details  our  presentation  so  far  follows 
that  in  [1]  and  [2]  •   In  this  section  we  supply  the  vali- 
dity proof  for  the  procedure,  which  is  missing  in  those 
accounts. 

We  shall  use  the  method  of  dominating  series:  First  a 
dominating  problem  is  constructed  whose  formal  solution  is 
given  by  series  of  positive  terms  that  are  term  for  term 
larger  than  the  absolute  values  of  the  corresponding  terms 
in  the  series  solution  of  the  original  problem.  Then  the 
formal  solution  of  the  dominating  problem  is  shown  to  con- 
verge, which  implies  the  convergence  of  the  original  series. 

We  begin  the  description  of  the  dominating  problem  by 
constructing  a  "dominating  function"  for  the  quantity  h  of 
(2.10).   Let  k  >  0  be  chosen  so  small  that  the  functions 
q(x),  r(x)  are  holomorphic  in  |x|  <  1  +  k  ,   If  K  is  an 
upper  bound  for  |q(x)|  and  |r(x)|  in  that  circle,  then  the 
series 

^(x  -  4)2  +  K^  (X  -  4)^  +  ...  =  I^(x-4)2(i.  ^^   )-l 

dominates    the  series    in  powers   of  x  -   4  of  the   two   expressions 
in  brackets   in   formula    (2.10).      Hence   the    function 
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"hi^/^,^-,   i)=^   (C-c)^(l-^)-^'^-l^)    (t.l) 


dominates  h  in  the  following  sense:  Expand  (2,10)  in  powers 
of  X  -  E,,    so  that  h  becomes  an  infinite  sum  of  terms  of  the 
form 


If 


^jmn^^)^^  -  C)*'  u  xt  ;   j  =  2^  3,  ...;  m  +  n  =  1 


^    ,  '^    » 1  ^m^  n 
Ajmn  (^  -^)^  ^^' 


is  the  corresponding  term  in  the  analogous  expansion  of  h, 
then 

|A,    (4)1  <  A,       In  kl  <  1  .         (1^.2) 

Now  consider  the   formal   series 


/\       ^       ^^  ^2 

X  =  4  +  x,a  +  Xptt     + 


x' =  1  +  X, a  +  XpQ     +    ... 

^        y\       /\  ^2 

u  =  u     +u,a  +  Uptt     +    ... 

u'=  u'    +u'a  +  ula     +    . . . 
o        1  2 


(U.3) 


where  x'  and  u'  are  not  the  derivatives  of  x,  u,  but  inde- 
pendent symbols  whose  numerical  value  is  irrelevant  at  this 
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point.      If   these    series   are    inserted   into  h  and   the    series 

rearranged  with  respect   to  powers   of   a,    then   the   coefficient 

h     of  a      is  a   polynomial    in   the  x.,    xl,    u.   which  dominates 
S  S  tr      ^  JJJ 

the  polynomial  h    introduced  before  in  analogous  fashion. 

Next  we   observe   that   the   fiinction 

(x   -   €)(u'    -  u'     ) 
o 

has  the  property  that  insertion  of  the  series  (U«3)  produces 

a  series  in  a  whose  s   coefficient  is  )       x  u  ,  which 

ti+v=s    ^     ^ 
H.,v>0 

is  analogous  to  the  first  summation  in  (2,v)»   Similar  ar- 
guments hold  for  the  remaining  terms  of  (2,v)«   Summariz.ing 
we  get  the  following  result:  Let  the  constant  K  introduced 
above  be  chOHen  so  large  that  also  |q«(x)|  <  K  in  |^|  <  1  +  k 
and  consider  the  function 


A  A  A  A  /^-        A    A  /^,    ^.-,  _^      ^,     ^  /., 
R(x,x«,u,u';^;a)=(x-^)(u'-u' )+  auu'  +  K(x'-l)(u-u  ) 


ih^k) 


+K(x-4)(x'-l)ui    +1    (x-^)''(k-(x-4))    "-(u+xM    . 
If  X,    x«,    u,    u'    are   replaced  by  the    series    (I|.»3)    then   the 

A 

coefficient  R  of  a^  in  the  resulting  series  in  powers  of  a 
is  a  polynomial  in  the  x.,  x.,  u.,  u.  with  positive  coef- 
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ficients,  independent  of  ^  and  not  smaller  than  the  absolute 
values  of  the  corresponding  coefficients  of  R  ,  in  I4I  <  1, 
for  s  =  1 ,  2 ,  .  . . 

In  order  to  avoid  an  accumulation  of  complexities  in  the 

presentation,  we  give  the  remainder  of  the  argument  first 

(.2) 
for  the  special  case  that  R^  '  =   0  for  all  s  .   The  modi- 

fications  necessary  if  R^  '  =|=  0  ,  for  1  <  s  <  m,  m  >  1,  will 

be  indicated  later.   For  our  special  case  we  have  u(4)  =  u  (4). 

Let  the  constant  c  be  chosen  so  large  that  in  addition  to 

the  inequalities  (3.15),  (3«l6)  the  relation 

!%(?)!  1   °f   °»  l%(^)l  <°€°    ,0<£<1     (i+.5) 


(cf,  formula  (3.U))  holds.   If  we  define,  furthermore,  u  by 


A        -^o 


we  can  take  as  dominating  problem 

A  ^o"^-^  ^  '^^   '^  -1   ^    '^   -1 

X  =  4  +  c^  ""   R  (X,  xc  ,  Uq,  Uq£  "■;    C,  a)   .      (1|.6) 


A  A 

In  fact,  if  X  is  replaced  by  £  +  ax,  +  ...  in  (1+.6)  we  find 

^«  =  °4  °   R.    » 
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where  in  the  polynomials  R   the  quantities  x'.  must  be  re- 

.  ^     -^  ^-1  ""  ,    ^"'^o  ''\    ^        J'^o~'^   and  u  ,  u '  ( j  >  0) 
placed  by  x^C   ,  %  by  c^    ,  u^  by  c?  J    j 

by  zero.   In  view  of  the  dominating  property  of  H  and  of  the 

inequalities  (3.l5),  (3.l6)  and  (U.5)  it  follows  inductively, 
that 

I  X  (f)|  1  X  .  0  <  4  <  1  .  (1+.7) 


Thus  the  convergence  of  the  series  for  >:(£)  will  be 
proved,  if  we  show  the  convergence  of  the  series  4  +  x-.a  +...  , 
or,  in  other  words,  if  we  can  show  that  equation  (U.6)  de- 
fines X  as  a  holomorphic  function  of  a  near  0,  =  0.   Written 
explicitly  (i+.b)   reads 


^  _  f  =  cc  °   [ac^^   °   +  K(x  -  4)(x  -  4)c£   ° 

+  |(x  -  4)^(k-(^-c))  ^c^  °+x4^)], 


or,  if  we  set  x-4=z,   a.^    =^» 


.3..  .2.JE    K^lc  ,„_  J«o.  ,Vl 


/\ 


c-^^  +  C%2  '^  +  ^  £_^  (c^  +  Z4  °  +  4  "   )  -  z  =  0  .    (U.8) 
The  left  aide  vanishes  for  z  =  0,  ?  =  0  .  At  this  point 

A 

its   partial  derivative  vjith  respect   to   z    is    different   from 
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zero.  In  fact  it  is  equal  to  -1  .   Hence  (U.S)  can  be  solved 
for  z  leading  to  a  holomorphic  function  of  ^  =  a^     near 
^  =  0  with  a  radius  of  convergence  that  is  bounded  away  from 
zero,  for  4inO<4<l*  This  proves  that  the  series  for 
x(4)  converges  for 


where  y  is  a  positive  constant  independent  of  4  «  VJe  suj:nnia- 
rize  our  results  so  far  as  a  theorem. 
Theorem  1,   In  the  differential  equation 


(x  +  au)  ^  +  q(x)u  -  r(x)   =  0 


the  functions  q(x)  and  r (x)  a r e  assume d  to  be  holomor phic  for 
lx|  <  1  .   Let  u(x)  be  the  solution  for  which  u(l)  =  b  .   If 


Qq  =  q(o)  >  0  , 


and  if  the  solution  u  (x)  of  the  problem 


du 
satisfies  the  condition 
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q(x)   u^(x)    -   r(x)   ^   0        in   0<>:<1 


then 


u(x)  =  u^(^)   ,  (U.IO) 


where  x  and  g  are  related  by  the  equation 


2 
X  =  f  +  x-,(f;)a  +  Xp(i^)c  +  ...  ,  (ii,  11) 


whose  coefficients  x  (^)  can  be  calculated  by  formal  in- 
sertion  of  (Ii-.lO)  and  (ij.ll)  into  (i|.9)  .   The  series  (U.ll) 
converges  for 


q 
|a|  <  Y^  °  ,    0  <  C  1  1  (1^.12) 

where  y  is  a  certain  positive  constant. 

Now  we  consider  the  more  t^eneral  case  that 

u(?)  =  u^(€)  +  u^(f)a  +  ...  +  u^(?)ct"'  (U.13) 


where  the  u.(£)  have  been  determined  by  solving  {2,1k)   with 
some  admissible  choice  of  Rj  '  .   This  means  that  in  our  ge- 
neral procedure  R^  '^  =  R   for  s  >  m, 
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The  constant  c  can  be  chosen  so  large  that,  in  view  of 
(3.11),  (2.13)  and  (2.1k) 


-jq  -(j+i)q 

|Xj|  <  c  kl    °,  lu.|  <  ck,|       °,  Id  <  1,  i    <   m,  (U.IU) 


and  also 

ix;i  <  cuT^'^^"  ju'.i  <  cui"  ^^  '^^^  ,ldfi,  J  <  m,  (U.15) 


Hence,  if  we  define 


-J^^o   A      -(0+1)^0 
Xj  =  c^     ,  Uj  =  C4         ,     j  =  l,...,m  ,      (i4-.16) 


we  have 


|Xj|  <  Xj  ,   |Ujl  1  Uj  ,   j  =  l,...,m;   0  <  C  <  1  •   (U.17) 


A 

Let  u  be  defined  by 


u  =  \    u  0.*^  ; 
0  =  0    j 

then 

A     -q„  T    ,   "^o.m+1 

u  =  c^  °  2^Z_i5^J .  (4.16) 
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/\        A    -1      ^        ^    -1 

Tith   the   definitions  x»    =  x4      ,    u'    =  uf      ,      i.e., 


x'    =  X  4"^    ,    u!  =  u  r   -^   ,        s   =  1,2,...,  (i|.19) 


it    is    furthermore   true    that 


|x!|    <   x!    ,        |u    I    <   u.'    ,  j    =    l,...,m;    0   <   £    <    1.     (i|.20) 

J  J  J  J 


Setting,  as  before, 

z  =  X  -  £,     ^  =  a4     , 


^»   /•  „     ""^o 


the  function  R  of  {k*k)  becomes  now,  after  some  simplifi- 
cations, 

A.  -q-l(  TviTi  o/i    ^x(\+l\2 

^   ,    ^m+1        ^  ';}2     /       ,    >,m+l        ^  q^  "^r."*"^';  t 


A 

However,  this  is  not  a  correct  choice  of  R  to  use  in 

(^.6)  in  the  present  case,  since  the  values  obtained  for 

X-, X  in  this  manner  would  be  different  from  those  de- 

1'    *  m 

fined  in  (L|..l6)  and  might  fail  to  satisfy  the  first  inequa- 
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lity  in  (Ij-.l?)  .   In  order  to  construct  a  dominating  function 
with  the  desired  properties  assume  first  that  m  =  1  and 
let  <1)-|(?)   be  the  function  obtained  for  x-,   instead  of 
c^     from  the  equation 

z  =  c?  °    R,  ('+.22) 


A 

R  being  defined  by  (a. 21) .   Then  the  same  calculation  re- 


peated with 


z  =  0?  °    R  +  a(cg   °  -  ^^{E,))  (1^.23) 


-^o 


does  yield  for  x,   the  required  value  OE,        ,      The  f miction 
<J)-,(?)  is  of  the  form 

^^i-)   =  ^"%'i"(.-)  (I;.2U) 

with  <|>n'(^)  bounded  in  0  <  ^  <  1.   This  can  be  shown  by  the 

A 

same  argument  as  that  leading  to  Theorem  1,  with  R  now 
defined  by  (I|..21)  .   Hence,  equation  {h.,23)    can  be  vjritten 

..  z  =   cs'^^'^^R  +  (c-f^-(§))^.  (1^.25) 
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-2q 
The  value  9p(€)   obtained  instead  of  c^      for  x-  from 

-2q  ,, 
this  equation  is  of  the  form  (|'p(?)  =  E,  <t'p(5)   with  bounded 

4p('?)  •   I^  order  to  obtain  the  correct  value  for  Xp,  formula 

(I4-.25)  must  therefore  be  replaced  by 

z  =  c^  °  R  +  (c-(J.i(?))^  +  (c-(t)|(c))^^.       (i!-.26) 

This  change  does  not  affect  the  value  of  x,.   Continuing  in 
this  fashion  an  eauation 


z  =  c^  °  R  +  EI  (c-<l>T(5))^^ 


is  finally  obtained  that  is  a  dominating  equation  for  our 
problem.   This  proves 

Theorem  2.   Theorem  1  can  be  generalized  as  follows; 
The  solution  u ( x )   can  be  represented  parametrically  in  the 
form 


u(?)  =  u^(5)  +  u^(^)a  +  ...  +  u^(.^)a"^ 


x(5)  =  ^  +  y:^{t,)a  +  yi^S^)^^   +  ..., 


vjhere  m  is  an  arbitrary  non-negative  intef^er.   The  coeffi- 
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clents   u  (•^),  X  (.^)   are  to  be  calculated  according  to  the 
formal  procedure  explained  in  sections  2  and  3,  setting 
R^2)  _  Q  ^^  s  >  m. 


The  series  for  x(^)   converges  for 


^o 


Y       being  a  certain  constant, 


al  <  Ym?  »        0  <  ?  <  1 


Corollary;   A  slightly  weaker  variant  of  the  statement 
of  Theorem  2  consists  in  saying  that  any  adinissable  vray  of 
splitting  R^   into  R;  '   and  R;  '   in  applying  the  method 
of  sections  2  and  3  leads  to  series  that  are  asymptotic  to 
the  solution  ^{£,) ,   u ( ^ )   of  the  differential  equation  In  the 
sense  that 


s=l  ^  ^ 


s=o 


with  E  (2, a)   and  P  (£,a)   uniformly  boimded  for 
———       m  ^ '    — —   m  ^ '    — — — ' 
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5.   Supplementary  Remarks 

a)   The  point  x  =  0.   The  main  purpose  of  Lighthill's 
construction  is  to  find  an  expansion  that  can  be  used  in  the 
neighborhood  of  x  =  0,  a  =  0.   In  order  to  show  that  this  is 
the  case  i-re  return  to  the  form  in  xrhich  our  main  result  was 
first  obtained.  We  had  found  a  function 

z  =  F(?,^)   ,  (5.1) 

where   z  =  x  -  ^,  holoinorphic  in  ^  for   1^1  <Yv,,j  O^C^l* 

The  number  y  i^  independent  of  %,      The  dependence  on  ^ 

is  holomorphic  in  0  <  E,   <  1     and  continuous  in  0  <  ?  <  !• 

Equation  (5«1)  is  satisfied  for   z  =  ^  =  0.   Hence,  the  in- 
verse function 

K   =  G(5,z)  (5.2) 

exists  and  is  a  holomorphic  function  of   z   in  a  neighborhood 
of  z  =  0,  for  values  of  .%      in  0  <  ?  <  1  such  that 
[BP(£;,^)/9^]^_Q  =  x'^'(:^)  7^  0.   If  complications  are  to  be 
avoided  we  must  introduce  the  hypothesis 


x'^(O)  ^   0  .  {S.3) 


si'        '.!. 


'••■'■•aC    ill 


-J-ij'x  ■-' 


•.'.:■  X.  '  :--r.'.   .     :<.- 


.  o  or  f  r- 


1  -.  ♦  ■ 


no  s./n; 


1    .-■  •  ■       ; 


iDOO.r 


;<-r 


-  25  - 

(The  case  x£'(0)  =0   is  briefly  discussed  in  [1],  p.  Il86)  . 
Then  xt{S,)    j^   0  for  0  <  ^  <  ?-j_  <  1.   Since  the  branchpoints 
^  vihere  3f/6^  =  0  depend  continuously  on  E,,    there  is  a 
number   y'   such  that  the  circle   |^|  <  yI*  contains  no  such 
branchpoint  for  0   <  ^  <  E,-,*     The  image  of  this  circle  depends 
continuously  on  E,     and  contains  the  point  z  =  0  in  its 
interior.   It  folloi/s  that,  although  the  circle  of  regularity 
of  the  function  G(|,z)   about  z  =  0  may  depend  on  E,,    it  has 
a  positive  lovjer  bound  5_,  for  0  <  E   <  E-,  .      The  series  for 
GC^^jZ)   converges,  in  particular,  for  z  =  -E,,   provided 

0  <  5  <  min(^3_,5^).  {S.k) 


If  ^   is  the  value  of  t^      in  (5.2)  corresponding  to  z  =  -E, 
we  get,  by  inversion. 


CD 


°      3=1   ^      ° 


i.e.,  X  =  0.   In  other  words,   x  vanishes  for  all  combinations 
of  values  of  E,     and   a   such  that 


r  '.-ii  ■ ',  r  f 


2Jaic>rir'orfo-^'' 


i     - 


4,..         ,     ^      .    .-.  i)jj 


J    f-!*"'         -i 


r*  *^'   *"   "'  ' 


'.     ■.-      ;    !  \ 


rrj  r   IC  ■;'.'■ 


K  r^   ,        r 


,  !  i  ;5  'A 


ilO 


(J.' ) 


ncl' 


c  I.  •■'■'. 


-(.ir  « 


-26'  -■ 


and  that  (5«U)  is  satisfied. 


,  00 

1  .    ■<: ,  ^  N_S 


Since  G(4,z)  =  — z  +  zH  8^(4)2  (g^(?)  holomorphic 

m  4  and  £;   ), 

we  have 

CO 

s=2 

and  therefore,  if  a  and  4  are  related  by  a^    ~  ^o  * 

,  q   cx) 

^  s=2 


The  right  member  is  a  once  continuously  diff erentiable  function 

^  i.1-     :.  X--1   4-    ^'^o"''   near  t  =  0.   By  an  application  of 
of  the  variable  t  =  ^  j  t^t- 

the  implicit  function  theorem  we  conclude  that 

^  ^   =-  x£(o)a  +  0(a) 

i.e* 

1/q^+l ) 
C  =  (-  x'£(o)a)    °    +  0(a).  (5.6) 


This  value  of  4  can  be  taken  real  if  x'£'(0)  <  0,  and  it  fol- 
lows that  our  series  for  x(^)  does  indeed  assume  the  value  zero 
for  a  certain  positive  £,,    if  a  is  small  enough.   If  x'£'(0)  >  0, 
Jhe  solution  x(4)  of  our  differential  equation  problem  does 
not  assume  the  value  zero  in  0  <  4  <  1.   Further  discussion 
of  this  situation  can  be  found  in  [1]. 
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b)   The  initial  condition.   In  an  example  given  in  [1] 
and  [2]  our  rigid  choice  of  the  boundary  conditions  as 

Xg(l)  =  0,   Uq(1)  =  b;   Ug(l)  =  0  ,     s  =  1,2,..., 
is  replaced  by  the  more  general  one: 


00         s  °°  s 

Cn  +  ^    ^s^^l^^  "  ■'■»    HI  ^g(^l)^  "  ^' 
s=l  s=o 


which  leads  to  a  recursive  set  of  relations  between  the  x^(^^) 
and  u  (E-,)*      This  method  has  the  advantage  that  half  of  the 
initial  values  can  be  chosen  arbitrarily  in  the  course  of  the 
calculation  so  as  to  simplify  the  resulting  expression.  An 
extension  of  the  argument  in  section  Ij.  which  would  include 
this  possibility  can  very  likely  be  given,  but  it  will  not  be 
attempted  hore. 

c)  A  limitation  of  the  method.  It  is  stated  in  [1]  that 
the  method  extends  without  essential  modifications  to  diffe- 
rential equations  of  the  form 

[x  +  ap^(x,u)  +  ...]  ^  +  q(x)u  =  r(x)  +  ar^(x,u)  +  ...  . 

I  fail  to  see  how  t..^:.fi  Ou,n  be  maintained  in  such  generality. 
If,  for  instance 
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r^(x,u)   =  u^ 

one  finds  -.„ 

%  "  %%  +  u^     =  0(4  )    , 

-2q   -1 
instead  of  0(4  )•      Thus  a  non-trivial  modfication  of   the 

procedure   seems   to  be  necessary.      Lighthill's   claim   is   cor- 
rect       if  p    (x,u)   and  r    (x,u)   are  polynomials    of  degree   s   at 
s  s 

most  in  u,  but  this  would  exclude  an  equation  as  simple  as 


/       V  du 

(x  +  au)  —  +u=asinu. 

dx 


6.   The  Convergence  Proof  for  q  <  0.  Leaving  aside  the 
case  q  =  0  in  which  logarithmic  terms  occur  in  the  series 
expansion  but  which  offers  no  essential  new  difficulties 
(cf,  [1]  and  [2]  we  assume  now  that 


q^  =  -  <<  0  . 
o 


The  convergence  of  Lighthill's  procedure  in  this  case  could 
very  likely  be  proved  by  the  methods  of  this  article.   Instead 
we  shall  treat  this  problem  by  a  variant  of  Lighthill's  tech- 
nique that  avoids  some  of  the  difficulties  he  encovmters, 
Fvery  solution  of  the  system 
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^x  =  jc  +  ay 

(6.1) 
^u  =  r(x)  -  q(x)u 


satisfies  the  single  equation  (1.1).   The  converse  Is  not  true, 
By  solving,  as  we  small,  this  system  instead  of  (1.1)  we  lose, 
therefore,  some  of  the  flexibility  of  Lighthill's  method. 
This  flexibility  could  be  regained  by  multiplying  the  right 
members  of  (6,1)  by  some  function 

P(t,a)  =  1  +  F-j_(t)a  +  P2(t)a  +  ... 

and  by  determining  the  P_(t)  successively  in  the  course  of  the 
calculation.   (Compare  the  appendix  by  Frdelyi  in  [2]  in  this 
regard).   But  we  shall  treat  only  equation  (6.1)  as  given. 

As  pointed  out  by  Lighthill,  there  is  some  simplification 
of  the  argument,  if  the  property 

r(0)  =  0  (6.2) 

is  secured  by  means  of  a  preliminary  translation  of  the  vari- 
ables.  This  will  be  assumed  to  have  been  done. 

If  ^  is  not  an  integer,  the  solution  u^(^;)  of  the  dif- 
ferential equation  problem 

4%  =  r(4)  -  q(^)u^,      u^(l)  =  b 
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is  then  of  the  form 


Uo(^)  =  4^Pi(^)  +  CP2(€)  =  4V"(C)  (^-3) 


where  p-,  (C)»  Vpi^)   are  holomorphic  for  0  <  ^  <  1, 

p  =  mini  K,    1),  (6.U) 

and  p'''(C)  holomorphic  in  0  <  C  5  1>  continuous  in  0  <  C  ^  !• 
If  /C  is  an  integer  a   term  of  the  form  E,     log  E,   appears. 
The  proof  is  not  seriously  complicated  by  this  term.   In  the 
interest  of  simplicity  we  shall,  however,  assume  that  t\ 
is  not  an  integer. 


Setting 


u  =  Uq  +  y  '■      (6.5) 


the   system    (6.1)    is   converted   into 

^x'    =  X  +   <i%(^)   +   ay  (6.6) 

^y'  =  -q(?)y  +  R(x,y,5j  (6.7) 


with 


R(x,y,C)  =  r(x)   -r(c)  -   [q(x)   -  q(4)][u^  +  y].    (6.8) 


We    insert 
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into  (6,6)  and  (6,7).   This  leads,  in  the  usual  manner,  to 
the  recursion  formulas 

^Xg  =  Xg  +  Ug_3_,    x^(l)  =  0,     s  >  0  (6.10) 

4Ug  =  -qu  +  1^,     u  (1)  =  0,     s  >  0  (6.11) 


where 

^s  ^  ^s(^l'*"»^s'^l'"**^s-l'^^  •  (^'l^) 

The  solutiors  of  (6.10)  and  (6.11)  are 

^s(^)  =  «/i  r  "S-l^^)  dt  ,  (6.13) 

^3(4)  =  C'^P(4)[  t"^'  ^p"*^(v)Vt)  dt.  .   (6. Ill) 


Here  p(£)  is  holomorphic  and  different  from  zero  for  0  <  4  <  1 

and  h  (f )  is  an  abbreviated  notation  for  the  value  of  h  in 
s  ^  "^  s 

(6,12)  at  f  =  '^  .    Hence,  we  have  the  appraisals 


1^3(4)!  <  c4^^^^1^^^|h^(^)I,   s>  0,   0<4<1    (6.16) 


with  certain  constant  c,  which  we  choose  so  large  that,  in 
addition,  (6.3)  implies 
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I%(C)I  <  C4P  .  (6-17) 

The  last  three  Inequalities  form  the  basis  of  a  conver- 
gence proof  for  the  series  (6,9).   The  fvinction 

R(x,y,4)  =  K  ^_\^^^)    (1  +  C4P  +  y)  (6.I8) 

is  readily  seen  to  dominate  the  function  R(x,y,4)  of  (6,8). 
Therefore  the  equations 

X  =  ^  +  a(ccP  +  9)  (6.19) 

u  =  c^  R(x,y,c)  (6.20) 

constitute  a  problem  that  dominates  the  differential  equation 
problem  (6.6),  (6.7)  withx(l)  =  1,  y(l)  =  0  in  the  previously 
explained  sense.   Insertion  of  (6.19)  into  (6,2o)  and  use  of 
the  definition  (6.I8)  yield 

u  =  c^^   a(c€P+u)     (i  +  ccP-Hu).  (6.21) 

k  -  a(c5^-^+u) 

This   equation  defines   u  implicitly  as   a  holomorphic   function 
of  the   variable 

'Yl  =   a^P  {6.22} 


n'x  .  ;^ii..i  .    S' 


J ;..    e-'.:vi  i  .  s>.fp;-i,  1    .:->'!. 


';ri!''^    ai'y    3;L.:i;''-;-:c'     ,'.f   n. 


<.Jj    :"' )■ ,   ■"'  ^    ■■■r.  )    '■'■ 


;  .'    i- 


,V!'-::"j''-.''y    o-u-     ■•s:;,.£;J:':!o''5    :r.-.ur   Ki 


)■      --•; 


i^- 


I.  '  s  '   ; 


-t  -      V  ^.1-.;  «  ■-  ,'      .-  ..  u 


;.  ■  >  ;;t 


,-i  J'     '"•■     :'-' 


'■.'TiW''::'.  .i.CV'."',    A    ft  (,:     ,' i  J  ';  S '!.J..,'~Ii;J     ?■     ",'■.'-■■•".%_*.■'"'•■<■ 


-  33  - 
in  a  circle 

where  5  can  be  taken  independent  of  ^.  The  relation  (6.I9) 
shows  then  that  x  -  ^,  too,  is  holomorphic  in  v?.  Thus  the 
following  theorem  has  been  proved. 

Theorem  3-   If  q(0)   is  negative  and  r(0)  =  0   the 

solution  of  (x  +  a,^)d5E  "*"  'l(^)^  ~  ^(^)  ~  ^  ^^^ ^^  ^(•'■)  ~  ^  ^^^ 

be  represented  by  tvjo  series  of  the  form 

X  =  ?  +^_  X  (?)a^,   u  =  u  (5)  +  ^  u  (?)a^ 
s=l  ^  °      s=l 

whose  coefficients  are  calculated  by  formal  insertion  into 
(6,1) .   The  series  converge  for 

|a|  <  5^"P  ,   p  =  min(K:,l)  ,   0  <  5  <  1  , 

where   5  is  a  certain  positive  constant. 

Thanks  to  assumption  (6,2)  we  have  x„(0)  =  0.  Hence 
x(0)  =  0,  and  the  preceding  theorem  gives  us  a  parametric 
representation  of  the  solution  curve  valid  between  x  =  1,  u  =  b 
and  X  =  0,  u  =  0. 
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